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ABSTRACT 

In many symbiotic stars there appear broad wings around Ha, of which the formation 
mechanisms proposed thus far include a fast outflow, motion of the inner accretion 
disc, electron scattering and Raman scattering of Ly/3. We adopt a Monte Carlo tech- 
nique to simulate the Raman scattering of ultraviolet photons that are converted into 
optical photons around Ha, forming broad wings, and compute its polarization. Not- 
ing that many symbiotic stars exhibit a bipolar nebular morphology and polarization 
flip in the red-wing part of the Raman scattered O vi features, we assume that the 
neutral scattering region is composed of the two components. The first component 
is a static cylindrical shell with finite thickness; and the second component is a fi- 
nite slab that is moving away with velocity Vp — 100 km s^^ along the symmetry 
axis of the first component. The cylindrical shell component yields polarization in the 
direction parallel to the cylinder axis. The strongest polarization is obtained in the 
limit where the height of the cylinder approaches zero and the scattering region effec- 
tively becomes a circular ring. As the height of the cylinder increases, the resultant 
polarization decreases and becomes negligible in the limit of the infinite cylinder. The 
polarization near the line-centre is weaker than in the far wing regions because of the 
large Rayleigh scattering numbers due to the large scattering cross sections near the 
line centre. The receding polar scattering component produces strong polarization in 
the direction perpendicular to the cylinder axis. In the presence of a Ly/3 emission-line 
component with an equivalent width ^ 0.5 A, the polarized flux exhibits a local max- 
ima at A = 6578 A that corresponds to the receding velocity with GAvp relative to Ha. 
When the both scattering components co-exist, the polarization is characterized by 
weak parallel polarization near the line-centre and strong perpendicular polarization 
in the red part. We discuss the observational implications of our computation. 

Key words: line: formation — line: profiles — polarization — radiative transfer — 
scattering — binaries: symbiotic 



1 INTRODUCTION 

Symbiotic stars are generally known to be interacting binaries of a mass losing giant and a hot white dwarf surrounded by 
an ionized nebula which is responsible for various prominent emission lines (Kenyon 1986). In the spectra of symbiotic stars, 
there exist very unique emission features around 6830 A and 7088 A. Schmid (1989) proposed that they are formed via Raman 
scattering by atomic hydrogen of O vi 1032, 1038 doublet. When an O vi line photon that is more energetic than Lya is 
incident upon a hydrogen atom in its ground Is state, it subsequently re-emits an outgoing photon with de-excitation into 
either Is state or 2s state. In the first case, we have a Rayleigh scattering process with no change in the frequency in the rest 
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frame of the scatterer. The second case corresponds to a Raman scattering process, from which we obtain 6830 and 7088 hne 
photons from incident O vi 1032 and 1038 photons, respectively. 

The scattering cross-sections associated with these Raman processes are of order a ~ 10~^^ cm^ (e.g. schmid 198t ; 
Nussbaumer, Schmid fc Vogel 198£| ; Sadeghpour & Dalgarno 1992; Lee & Lee 1997), and therefore the operation of the O vi 



Raman scattering requires the existence of both a highly ionized O vi emission region and a fairly extensive neutral scattering 
region with H i column density Nhi ^ 10^^ cm~^. The scattering cross section is very sensitive to the wavelength of the 
incident radiation in such a way that it steeply increases as the wavelength approaches the resonance wavelength of the 
scatterer. Therefore, continuum photons with wavelength near that of the Ly/3 will be scattered in a neutral hydrogen region 
with much lower H i column density than that is required for O vi Raman scattering. 

van Winckel, Duerbeck & Schwarz (1993) showed that most symbiotic stars in the southern hemishpere exhibit very broad 
wings around Ha. A similar result was reported for those in the northern hemisphere by Ivison, Bode & Meaburn (1994). The 
full width at zero intensity often exceeds 2000 km s^^. Similar broad wings also appear in young planetary nebulae, including 
M2-9, Mz3, and IC 4997 (e.g. [Miranda, Torrelles fc Eiroa 1996|; [Arrieta fc Torres- Peimbert 2002|). In particular, for M2-9, 



Bahck (198!:) reported the existence of very broad Hq wings with a width ~ 10 km s . Remarkably, selvelli & Bonifacio 



(2000) presented the Very Large Telescope (VLT) spectrum of the symbiotic star RR Tel which shows Ha wings extending 
up to 12000 km s^^. According to van de Steene, Wood & van Hoof (200C), post-asymptotic giant branch (post-AGB) stars 
also exhibit broad Ha wings. 

There have been many theoretical suggestions for formation mechanism of the broad Ha wings. These include fast outflows 
(e.g. Schmutz et al. 1994), electron scattering (e.g. Bernat & Lambert 1975), the accretion disc (e.g. Quiroga et al. 2002) and 



the Raman scattering of Ly(3. In this paper, we will not consider all those theoretical models but will focus only on the Ha 
wing formation via the Raman scattering of Ly/3. 

Nussbaumer et al. (1989) discussed the astrophysical importance of Raman scattering by atomic hydrogen and proposed 
that broad Ha wings can be formed through the Raman scattering of Ly/3. This idea has been applied to young and compact 
planetary nebula IC 4997 by Lee & Hyung (2000). They computed the wing proflle formed through Raman scattering of flat 
continuum around Ly/3 and obtained an excellent agreement with the observed proflle. Similar success has been achieved for 
many symbiotic stars by Lee (2000). 

The existence of a thick H i region in symbiotic stars and young planetary nebulae can also be inferred from the He ii 
Raman scattered features. Because He ii is also a single electron atom, the wavelengths associated with the transitions 
from n = 2m state to the 2s state are very close to those corresponding to the mth Lyman series of hydrogen. Therefore, 
He II Raman scattering requires much less H i column density Nhi ~ lO'^" cm~^ than O vi 1032, 1038 doublet. The He ii 
Raman-scattered lines have been found in the spectra of RR Tel and He 2-106 (e.g. van Groningen 1993; Lee, Kang & Byun 
2001) and the young planetary nebulae M2-9 (Lee et al. 2001), which highlights the plausibility of the operation of the Ly/3 
Raman scattering in these objects. Therefore, in this work, we will consider Raman scattering of Ly/3 in neutral regions with 
Nhi = 10^'' - 10^^ cm'^ 

Spectropolarimetry will provide valuable information on the Ha wing formation mechanism. Raman-scattered radiation 
consists of purely scattered radiation without any mixing of direct radiation, and therefore Raman-scattered radiation can be 
strongly polarized. Polarization is also dependent on the scattering geometry. Spectropolarimetry may reveal the circumstellar 
matter distribution and be useful to understand the mass-loss process of the giant component. 

In this paper, we will compute the profile and the polarization of Ha wings that are assumed to be formed through 
Raman scattering of Ly/3. In section 2, we will briefly introduce the basic atomic physics relevant to the Ha wing formation 
and describe the Monte Carlo procedure. In section 3, we present our main results. In the flnal section we discuss our results 
and observational implications. 



2 Ha WING FORMATION THROUGH RAMAN SCATTERING 
2.1 Branching Ratio 

A Ly/3 photon incident on a scattering region consisting of atomic hydrogen undergoes either Rayleigh scattering or Raman 
scattering. In particular, in the latter case the scattering H atom de-excites to the 2s state resulting in an outgoing Ha photon. 
Owing to the scattering incoherency, if an incident photon near the Ly/3 line centre is Raman-scattered, then the outgoing 
photon finds itself in the Ha wing region much further from the line centre, by a factor of 6.4, than from the line centre of Ly/3 
(e.g. Schmid 1989, Nussbaumer et al. 1989). Lee & Hyung (2000) emphasized that the profile of Ha wings formed by Raman 
scattering of Ly/3 will be proportional to AA~^, when the Raman scattering optical depth is small. However, the detailed 
conversion process from a ultraviolet photon into an optical photon is dependent on the branching ratio of the Rayleigh 
scattering and Raman scattering cross sections (e.g. Schmid 1996). In this subsection, we illustrate a basic calculation that 
provides an approximate branching ratio near the Ly/3 line centre. 
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Figure 1. Branching ratio of the Raman scattering to the Rayleigh scattering. The solid and the long dashed lines represent the 
differential cross sections of Rayleigh and Raman scatterings in units of Thomson scattering cross section, respectively. The branching 
ratio of the Raman scattering to the Rayleigh scattering and its first-order approximation are shown in the dotted and the short dashed 
lines, respectively. 



The branching ratio of the Rayleigh process to the Raman process is weU known in the hterature, which is around 1:7. 
However, this ratio is dependent on the wavelength of the incident photon. 
The differential cross-section for Rayleigh scattering is given by 
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where the subscript A stands for the initial quantum state of the scatterer, / for the intermediate stat e and B for t he final 
state. We have A = ls,B = 2s for the Raman process and A = B = Is for the Rayleigh process ( Bakurai 1967 ). Here, 
^^lA ~ h~^{Ei — Ea) is the angular frequency for the transition between the initial state A and the intermediate state /, 
ro is the classical electron radius, and e" represents the polarization vector of the incident photon. The primed quantities 
correspond to those associated with the outgoing photon. In particular, the energy conservation relation for the Raman 
scattering case gives 

Uj' = U) — UJa, (3) 

where uic, represents the angular frequency of the Lya transition. We are interested only in scatterings occurring in the wing 
regime far from the resonance and therefore omit the damping term in the denominator of each transition matrix element. 

For an incident photon around Ly/3, u> ~ u)0 and the cross sections for both scattering processes are predominantly 
contributed by the first term in the summation corresponding to / = 3p in equations (jl|) and (P). Since the angular 
distribution of the scattered radiation is the same for both processes, we can write the branching ratio rj, of the cross-sections 
for Rayleigh scattering and Raman scattering as 
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where we set w — uj0 + Acu, \Aijj\ iop. 

The matrix elements that appear in equation 



are explicitly given by 



< 2s II p II 3p >< 3p II p II Is >, 
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Figure 2. Schematic view of the scattering geometry. The central core represents the emission region. The static cyhndrical shell denoted 
by Component I is composed of three layers stratified from the equatorial region. The layers are characterized by different H i densities 
and heights hi,h2, h^. The circular slabs (Component II) with homogeneous H i density are located in the polar direction, moving away 
from the core. 
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It is noted that the summation extends to the continuum states, where the sum over all continuum states should be written 
in an integral form. 

The explicit representations of the matrix elements are found in the literature (e.g. [Karzas fc Latter 1961; Saslow & 



Mills 1969 



Berestetskii, Lifshitz 



Pitaevskii 1971) and a straightforward numerical computation gives Ai/{ujpA3) — 0.9268, 

to yield 



-30.37. These results are directly substituted into equation 



ujf3A4/As = 17.91, uj^A2/Ai = 

rt{uj) ~ 0.1342 - 12. 50 Alo/uj p. (6) 

In Fig. ^ we present the branching ratio of the Raman scattering cross-section to the Rayleigh cross section. The short- 
dashed line represents the approximate result given in equation (^), whereas the dotted line stands for the result obtained by 
using the whole Kramers-Heisenberg relation. The approximate relation gives an error of less than 10 per cent in the range 
AV ^ 1000 km s~^. This implies that the contribution from other np states is not negligible. From Fig. |l| we see that the 
branching ratio varies in a rather large range from 0.05 to 0.45 in the wavelength interval from 1015 A to 1040 A. Since 
the broad feature around 7088 A which is Raman-scattered O vi 1038 is often found in the spectra of symbiotic stars, it is 
necessary to investigate the Raman conversion efficiency for a broad range of the branching ratios. 



2.2 Scattering Geometry 



Many symbiotic stars exhibit a bipolar nebular morphology, which is characterized by the accumulation of circumstellar 
materials near the equatorial plane. Many studies of Raman scattered O vi lines have demonstrated that the neutral scattering 



region can be associated with the cool component, while the ionized region must be close to the hot component (Harries & 



Howarth 1997; gchmid fc Schild 1997; gchmid, Corradi, Krautter fc Schild 2000) 



However, the Raman scattering of Ly(3 becomes significant in neutral regions with Nhi ^ lO^^cm"^, which is two orders 
of magnitude smaller than the column density required for the operation of O vi Raman scattering. This implies that the Ha 
wing formation through Ly/3 Raman scattering should occur in much more extended regions than those that can be probed by 
Raman-scattered O vi features. This is consistent with the fact that most symbiotic stars exhibit broad Hq wings whereas only 
about half of them show Raman-scattered O vi features (van Winckel et al. 1993, Ivison et al. 1994). Furthermore, broad Ha 
wings are also found in young planetary nebulae and post - AGS stars, which appear to hav e common characteristics including 
bipolar morphology and composite emission- line profiles ( Arrieta fc Torres-Peimbert 2002 ). It is still uncertain that all these 
objects are binary systems, whereas symbiotic stars are generally regarded as binary systems. In view of these circumstances, 
we adopt a simple neutral scattering region that is approximated by a cylindrical shell. 

Fig. ^ shows a schematic view of the scattering geometry. We just assume that the neutral hydrogen component forms a 



Polarization of the Broad Ha Wing in Symbiotic Stars 5 



finite cylindrical shell (Component I in Fig. with the centre coinciding with the Ly/3 emission source. The cylindrical shell 
is assumed to be of constant H i density and hence characterized by radius r and height h. 

In the spectropolarimetric data presented by Harries & Howarth (1996), the Raman-scattered O vi 6830 and 7088 features 
exhibit polarization flip in the red-wing part. 3chmid {199t ) interpreted the polarization flip in the red-wing part by scattering 
in the slow stellar wind from the giant component. Another interpretation has been proposed by Lee & Park (1999), who 
assumed the existence of an accretion disc around the white dwarf component via gravitational capture of the slow stellar 
wind around the giant component. According to their model, the emission region is formed around the accretion disc, and 
the main double peaks shown in the Raman-scattered features are attributed to the disc motion viewed in the equatorial 
direction, where the main scattering region is formed around the giant component. In order to explain the polarization flip in 
the red-wing part, there has to be another scattering component that is moving away from the emission region in the direction 
perpendicular to the equatorial plane. 

In this work, we will also consider that Ha wings are significantly contributed by this polar scattering region (Component II 
It is uncertain whether this polar scattering component is ubiquitous in all symbiotic stars. Harries & Howarth 



in Fig. I 
(1996[ 



Park 



showed that a significant fraction of symbiotic stars may possess this component, in the picture presented by 



Lee & 



199S:). Because the polar component is assumed to be receding from the emission region, it will mainly affect the red 



part of the Ha wings. 

With this geometry, we prepare a photon source at the centre, from which the ultraviolet photons around Ly/3 are 
injected into the scattering region. For a given incident photon the corresponding cross-sections of the Rayleigh and the 

|Bethe 



fc Salpeter 1957| ; [Sakurai 1967| ). The 
Lee 1997). According to this ratio, we 



Raman scatterings are calculated from the Kramers-Heisenberg formula (e.g. 
ratio of the two cross-sections is dependent only upon the incident wavelength ( Lee 
decide whether a given scattering event is Rayleigh or Raman in the Monte Carlo procedure. 

Although the ultraviolet photons may be Rayleigh-scattered several times, the scattering regions are assumed to be 
transparent to the optical photons converted by Raman scattering. Therefore an incident ultraviolet photon suffers a number 
of Rayleigh scatterings until it is Raman scattered or escapes the scattering region as an ultraviolet photon. 

We typically inject 10^ ultraviolet photons into a bin with AA = 0.005 A around Ly/3, and calculate the polarization 
state at each scattering. Finally we collect Raman scattered photons according to the emergent direction cosine of the wave 
vector with a bin size of A/i = 0.1. 



2.3 Computation of Polarization 

Because our scattering geometry is azimuthally symmetric, under the assumption that the emission source is not circularly 
polarized, no circular polarization will develop. Hence we only consider linear polarization that may develop either in parallel 
or perpendicular to the symmetry axis. We will distinguish the polarization direction by the sign of the polarization. That 
is, we will denote the polarization perpendicular to the symmetry axis by positive polarization and parallel polarization by 
negative polarization. 

In the case of single electron atoms including hydrogen, the fine structure level splitting is small because of accidental 
degeneracy. Therefore any off-resonance scattering is characterized by the same scattering phase function as that of the 



Thomson scattering, which was shown by Stenflo (198C). Since the Ha wings spread over ^ 1000 km s~ , resonance scattering 
contributes very little to the wing formation. Hence, in this work, the relevant scattering phase function is that for Thomson 
scattering. 

In the computation of linear polarization, we adopt the density matrix formalism, in which the polarization state is 
described by a 2 x 2 Hermitian matrix p. In particular, with respect to the symmetry axis which is taken as the z-axis in the 
usual spherical coordinate, we choose the polarization basis vectors by ei = (sin (j>, — cos (j>, 0), £2 = (cos 9 cos (j), cos 9 sin (j>, sin 9) 
for a given line photon with wavevector k = (sin S cos 0, sin 61 sin 0, cos 0), the linear polarization p is obtained from the 
difference of the main diagonal elements of the density matrix p, 

P = Pll - P22- (7) 

The relation between the density matrix p' of the scattered radiation with wavevector k = (sin 9' cos (j)' , sin 9' sin (j)' , cos 9') 
and p of the incident radiation is given by 



p'li — pii cos^ Aif) — pi2 COS 9 sin 2A(f> + P22 cos^ 9 sin^ A(j>, (8) 
p'i2 = ^Pii cos 9' sin 2A0 + pi2(cos 9 cos 9' cos 2 At/) + sin 9 sin 9' cos A^) — P22 cos 9{sm 9 sin 9' sin A(f> + i cos 9 cos 9' sin 2A(f>), 



p'22 = Pii cos"^ 9' sin^ A<jf> + pi2 cos 9' {2 sin 9 sin 9' sin A(f) + cos 9 cos 9' sin 2A(f)) + p22(cos 9 cos 9' cos A^ + sin 9 sin 9')^ , 
where A^ = (f>' — (j). 

It can be also shown that the angular distribution of the scattered radiation for an incident photon with k and p is given 



by the trace of p' (9' , (f)' , 9 , (f)) (e.g. Lee, Blandford fc Western 1994). Therefore, once the wavevector fe of a scattered photon 
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Figure 3. Polarization of the Raman scattered photons for various lines of sight. Ultraviolet photons are normally incident on an infinite 
slab with a finite thickness. Various types of lines represent the number of scatterings that outgoing photons suffer. The polarization is 
set to be positive when the polarization direction is perpendicular to the slab normal. The horizontal axis represents the cosine of the 
angle {/j. = cosd) between the slab normal and the line of sight. The error bars represent one standard deviation. 



is chosen from p' , the polarization is also determined at the same time. It is also noted that the density matrix is normalized 
each time so that it has a unit trace. 



3 RESULTS 

First, we present two idealized cases in order to obtain a qualitative understanding of the polarization development through 
Raman scattering. Next, adopting a simple geometry we calculate the profile and the polarization for various geometrical 
parameters. Then we present two applications for more realistic geometries. 



3.1 Infinite Slab 

To investigate the dependence of the polarization on the number of Rayleigh scattering that a photon suffers, we just focus 
on the scattering events, ignoring other effects. An incident photon can escape and be counted only if it attains the prescribed 
number of scattering. 

We typically prepare 10* photons which are normally incident on an slab with a finite thickness. Here, we assume that the 
horizontal physical dimensions of the slab are infinite. We collect only those photons that have been scattered the prescribed 
number of times. We present the result in Fig. where the horizontal axis represents the cosine of angle {p — cosO) between 
the slab normal and the line of sight. The polarization steeply decreases as the scattering number increases. If a photon is 
scattered three times, its maximal polarization reduces to about 40 per cent, and almost zero polarization is observed when 
the number of scatterings exceeds 10. The general behaviour of polarization weakening as a result of increased number of 



scatterings can be found in the literature (e.g. Angel 1969). 

Next, we present an idealized case of incidence normal to the slab in Fig. ^. We also inject 10* ultraviolet photons in a 
bin with AA = 0.05 A with full width 6 A centred at Ly/3. They are normally incident to the infinite slab lying in the x — y 
plane. The H i column density is set to be 5 x lO^^cm"^. 

Despite the large column density, for an incident ultraviolet photon far from the resonance wavelength of Ly/3, the 
scattering region is optically thin and hence the Raman conversion rate is very low. However, radiation in this wing region is 
highly polarized up to 100 per cent in the equatorial direction, because the scattering number rarely exceeds 2. 

As the wavelength goes close to the Ly/3 line centre, the cross-section increases steeply. In this case, most incident 
ultraviolet photons may suffer a number of Rayleigh scatterings before escape by Raman scattering. They escape from the 
slab by the final Raman scattering and contribute to the broad Ha wing formation. The flux near the Ly/3 line centre 
is typically polarized by about 40 per cent in the equatorial direction, which is weaker than in the far-wing regions. The 
weakening of polarization in this region is attributed to the increased Rayleigh scattering numbers before escape. 
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Figure 4. Profiles and polarization of Raman scattered photons for various lines of sight. 10* ultraviolet photons are injected normally 
on to the infinite slab. The H i column density is set to be 5 X lO^'^cm"^. The error bars represent one standard deviation. 



3.2 Polarization for Various Heights of the CyUndrical Scattering Regions 

We assume that the H i number density is constant throughout the scattering region of Component I in Fig. ^ (/i = /ii = /i2 = 
/la), and obtain the profile and polarization of broad Hq wings, changing the ratios hjr of height to radius. In this work, just 
for simplicity and with a lack of detailed information on the circumstellar matter distribution, we assume that the thickness 
of the H I cylindrical shell is much smaller than the radius. The profiles and polarization of the broad Ha wings for various 
lines of sight are shown for several ratios in Figs. ^ and ^ 

The radiation emergent in the polar direction is almost unpolarized because of the symmetry, which provides a check 
of our code. When the observer is in the equatorial direction, the Raman scattered radiation exhibits negative polarization. 
As the height-to-radius ratio approaches zero, the scattering region is regarded as a circular ring, which yields the strongest 
polarization, 33 per cent, in the direction parallel to the symmetry axis of the circular ring. Since photons falling on to far-wing 
regions are scattered fewer times than those contributing to near-wing regions, the polarization in far-wing regions is much 
stronger than that near the Ha centre. Therefore it is particularly notable that the Raman scattering wings are characterized 
by weak polarization near the centre and strong polarization at far-wing regions, in the direction of the symmetry axis. 

As the height increases relative to the radius, for a given line of sight, there can be a broad range of scattering angles. 
Therefore a superposition of electric vectors associated with the scattered radiation lessens the overall polarization. In Fig. ^ 
it is apparent that the polarization weakens as hjr increases. When /i/r = 2, we obtain Ha wings polarized up to 10 per cent 
observed from the equatorial direction, which is much smaller than 30 per cent obtained in the case of hjr = 0.5. 

This trend will continue until the geometry is approximately regarded as an infinite cylinder. In this case, when the cylinder 
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Figure 5. Profiles and polarization of broad Ha wings for various lines of sight. Here we consider only Component I in Fig. g with 
h = hi = h2 = hs. The height-to-radius ratios h/r o{ the H i cylindrical shell are 0.5 and 1 for left- and right-hand panels, respectively. 
10^ ultraviolet photons in a bin with AA = 0.005 A are injected to the H i component I. The homogeneous H i column density is set to 
be lO^'^cm"^. The error bars represent one standard deviation. 
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has enough thickness, emergent Raman-scattered wings will be contributed to equally from all directions of incidence, and will 
be completely unpolarized. We increase the ratio h/r to as high as 5, which corresponds to the case that ultraviolet radiation 
from the emission source with 9 > 10° will be incident on to the scattering region. As is shown in upper panels of Figs. ^ and 
we obtain negligible polarization in this case for all lines of sight. 
However, we also see the weakly polarized parts in far blue and red regions for the case h/r — 5. This polarization, 
developed in the direction perpendicular to the symmetry axis, is attributed to the limb effect. That is, the cylindrical shell 
is optically thin in the direction normal to the shell for far blue- or red-wing photons. Therefore, only a small fraction Tw of 
photons incident in this direction will be scattered, where is the Raman scattering optical depth in this direction. However, 
when the incident photon has the wavevector k with fc • i; = ~ 1, the fraction of scattered photons increases to T-wj \J\- — 
which is much larger than t™. This limb effect strengthens the positive polarization component and results in polarization in 
the perpendicular to the symmetry axis. 

In the cylinder scattering model, we may summarize by saying that the polarization develops in the direction parallel 
to the symmetry axis, and that the polarization is weaker near the Ha centre, because of large number of the Rayleigh 
scatterings, than it is in the far-wing regions. 



3.3 Hybrid Model of Cylindrical Shells 

We present a hybrid model of Component I in Fig |2| consisting of three stratified layers {h\ < hi < hz) with different H i 
column densities. In this model, the density is highest in the equatorial region and decreases as the latitude increases. This 
may mimic the matter distribution in the circumstellar region more realistically. This hybrid model can be viewed as a linear 
superposition which is the solid-angle-weighted sum of the previous models with various ratios and homogeneous hydrogen 
number density. 

We consider two different cases denoted by Model I and Model II. In Model I, we set h\ = O.lr, hi = 0.5r and h^ = 5r; 
in Model II, we set /ii = 0.5r, hi = r and /13 — 5r. In both cases, the H i column densities of the stratified layers are chosen 
to be Nhi = 10^\ 10^" and 10^'' cm'^ respectively. 

We obtain similar patterns in polarization, which is weak near the Ha centre and becomes strong in far-wing regions. 
These polarization patterns are also attributed to larger scattering numbers for photons near the line centre than those for 
far- wing photons. 

It is particularly notable that a flat profile near the line centre no longer appears. This is because of the contribution 
made by the small column density regions at high latitudes. These regions are optically thin for most of incident ultraviolet 
photons, and hence yields profiles proportional to the scattering cross-section which is very steep near the line centre. In 
Model II shown in the right-hand panel of Fig. ^, the profile is fairly extensive because of the much larger column density 
than in the case of Model I shown in the left-hand panel of the same figure. 



3.4 Raman Scattering from the Receding Polar Components 

In this subsection, we consider the Raman scattering from the receding polar scattering region of Component II, the existence 



of which is implied by Harries & Howarth (1996). In this work, we simply assume that this component is moving away from 
the emission source along the symmetry axis of the cylindrical shells considered in the previous subsections. 

In the left-hand panel of Fig. ^ we show the polarization of Ha wing photons Raman-scattered in the receding polar region. 
We set the recession velocity of the polar component to be Vp — 100 km s~^. Since the symmetry axis and the scattering plane 
make a very small angle, the polarization develops in the direction perpendicular to the symmetry axis, or in our formalism 
we obtain positive polarization. For an observer in the equatorial direction, the scattering angle is nearly 90°, and therefore 
almost complete polarization is obtained. 

The wing flux is proportional to the scattering cross-section and the incident UV flux in the optically thin limit. In the 
observer's rest frame, the scattering cross-section peaks at the recessional velocity of the polar component Vp, and therefore 
we obtain a small peak at this velocity. Hence the resultant profile is obtained effectively by shifting the wing flux that is 
formed from Raman scattering in a static medium. Owing to the incoherency of Raman scattering, the wing fiux peaks at 
6578 A. This is interesting in that the symbiotic stars RR Tel and He 2-106 appear to exhibit a small and broad bump around 



6580 A just blueward of [N II] 6584 feature ([Lee, Kang fc Byun 2001 ) 



In the right-hand panel of Fig. ^ we present the Raman scattered wing fiuxes obtained from combinations of a static 
cylindrical shell and a receding polar component. Since the two scattering components provide opposite senses of polarization, 
the resultant polarization is obtained by the algebraic sum of the two polarized flux components divided by the total flux. 
Therefore, in the blue wing, where the polar component contributes little, we obtain polarization in the direction parallel 
to the symmetry axis. However, in the red wing, the polar component predominantly contributes to the final polarization, 
yielding perpendicular polarization. 
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4 DISCUSSION AND SUMMARY 

We have investigated the polarization of Ha wings formed through Raman scattering of Ly/3. Two kinds of scattering regions 
are considered. The first one is a static cylindrical shell. From this component, the polarization of Ha wings develops in the 
direction parallel to the cylinder axis. It is particularly notable that stronger polarization is obtained in far-wing regions than 
near the line centre. Near the line centre, the scattering cross-section is larger, and therefore Rayleigh scattering occurs more 
frequently than in the far-wing regions. The weaker polarization in the near-wing regions is attributed to the larger Rayleigh 
scattering numbers than in the far-wing regions. 

We consider the second scattering component which is located in the polar regions and moving away from the ultraviolet 
emission source with a recessional velocity Vp = 100 km s~^. The physical nature of this component is not certain, but the 



existence of this component is apparent from the spectropolarimetric observations presented by Harries & Howarth (1996). 
This component produces polarization in the direction perpendicular to the symmetry axis. Since the scattering angle is 
almost 90°, nearly complete polarization is obtained. Since the scattering component is moving away from the source, the 
wing flux peak appears at ?; = 6Avp, or around 6578 A with our choice of Vp — 100 km s~^. Therefore, this polar component 
will produce a small bump in the wing profile which is strongly polarized in the direction perpendicular to the symmetry axis. 

When the both scattering components co-exist, we may observe the polarization flip, where the blue part is characterized 
by parallel polarization and the red part is polarized in the perpendicular direction. Therefore, if the scattering region is 
composed of a cylindrical shell subtending a fairly large solid angle with sufficient column density and a small receding 
polar component, the overall wing profile is determined by the cylinder component and the polarization will be predominantly 
contributed by the polar component in the red part. Therefore in this case, we mainly obtain a profile that is nearly proportional 
to AA~^ and strongly polarized in the red part in the direction perpendicular to the symmetry axis. 



This is very interesting, considering the spectropolarimetric observation of the symbiotic star BI Crucis by Harries (1996), 
who reported the polarized Ha wings in the direction perpendicular to the symmetry axis. They also noted that the polarization 
is very strong in the red part. In their spectropolarimetric data, no polarization flip is apparent. They interpreted their data 
using the assumption that the wings are formed through electron scattering. However, their data appear to be qualitatively 
consistent with the hypothesis that the Ha wings are formed by Raman scattering in a fairly extensive cylindrical shell and 
a receding polar component. 

Lee (1999) proposed that electron scattering of monochromatic line radiation may leave polarization structures. The 
radiative transfer process through Thomson scattering is diffusive both in real space and frequency space. This implies that 
photons that suffer more electron scattering tend to fall in both far- and near-wing regions, whereas those with fewer numbers 
of scatterings contribute to the formation of near-wing regions. Since polarization becomes weaker as the scattering number 
increases, we obtain weaker polarization in the far-wing regions in Thomson scattering media. Therefore spectropolarimetry 
can be an important tool to investigate the wing formation mechanism. 

The Ly/3 photon source was considered to be point-like. In reality, this photon source will be extended and even partly 
coincide with the scattering region. This will most likely lower the resulting polarization. The Ha line flux will be strongly 
diluted by photons from the emission nebula, originating predominantly from the recombination process. This contribution 
can be dominant in the line centre and it will lower the resulting polarization considerably. Other factors that may affect the 
polarization of Raman-scattered Ha wings and have not been considered in this work may include dust and the kinematics 
of the cylindrical shell component. Currently, only a few spectropolarimetric observations of Ha wings can be found in the 
literature, and little is known about the formation of Ha proflles seen in symbiotic stars. We hope that more observations will 
provide constraints on various theoretical models of Ha and other Balmer series line formation. 
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Figure 6. Profiles and polarization of broad Ha wings for various lines of sight. The height-to-radius ratios of the H i cylinder are 2 and 
5 for left- and right-hand panels, respectively. Other parameters are the same as described in Fig. 
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Figure 7. Profiles and polarization of broad Ha wings for the liybrid model. We consider only Component I in Fig. ^. The H i column 
densities are set to be 10^^, 10^" and 10^^ cm~^ from the equatorial region with hi, h2, /i3=0.1, 0.5 and 5 in units of the radius r for the 
Hybrid Model I, and 0.5, 1 and 5 for the Hybrid Model II. Other parameters are the same as described in Fig. pi 
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Figure 8. Profiles and polarization of broad Ho wings for both components. We set the recession velocity of the scattering region of 
Component II to be Vp = 100 km s~^, the half-opening angle to be 20° and the ratio of Component I (h = hi = h2 = hs) to be h/r = 0.5. 
10® ultraviolet photons in a bin with AA = 0.005 A are injected into only Component II with H i column density 10^" cm~^ in the 
left-hand panel, and into both components (H i column density 10^", 10'^'' cm~^ for Component I and Component II, respectively) in 
the right-hand panel. We omit to denote the polarization when the photon flux is negligibly small. The error bars represent one standard 
deviation. 



